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ABSTRACT 

The  mcdel  presented  in  this  thesis  is  a  computer  simulation  model 
of  ASW  interactions  between  a  formation  of  high,  value  group  ships, 
protected  by  some  screening  ships,    and  some  penetrating  submarines. 
The  model  is  designed  for  use  as  an  aid  in  improving  the  ability  of  a 
proposed  screening  tactic  in  the  detection  of  a  penetrating  submarine. 
A  systematic  procedure  to  improve  a  screen's  effectiveness  against  a 
known  submarine  threat  is  demonstrated,    and  an  example  problem  is 
worked  using  this  procedure. 
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I.     INTRODUCTION 

In  the  present  era  of  a  continuing  and  growing  Sino-Soviet  sub- 
marine threat,    the  problem  of  effective  anti-submarine  warfare  (ASW) 
to  protect  a  high  value  group  of  shipping  is  of  considerable  interest  to 
the  Naval  Officer.     In  recent  years  a  great  deal  of  effort  has  been 
expended  in  an  attempt  to  develop  an  optimal  solution  to  the  ASW  prob- 
lem; but,    due  primarily  to  the  complexities  of  the  environment  and  the 
resulting  inability  to  quantitize  parameters,,  no  such  optimal  solution 
has  been  achieved.     It  is  not  the  purpose  of  this  paper  to  provide 
another  solution  technique  for  the  ASW  problem,    but  rather  to  provide 
a  computer  simulation  model  of  an  ASW  action  and  to  suggest  one  pos- 
sible use  for  the  model. 

The  model  presented  in  this  paper  is  designed  for  use  as  an  aid  in 
the  design  of  ASW  screens  by  providing  measures  of  effectiveness  of  a 
screen  in  the  detection  of  penetrating  submarines,    the  number  of  attacks 
made  by  penetrating  submarines,   and  the  number  of  successful  attacks 
made  by  a  penetrating  submarine.     Using  the  computer  simulation 
model,    the  effectiveness  of  a  simple  ASW  screen  against  penetrating 
submarines  is  maximized  within  the  constraint  limitations  provided. 


H.     DESCRIPTION 

The  model  is  a  computer  simulation  of  the  ASW  environment, 
utilizing  an  event  store  procedure.      The  model  is  written  in  FORTRAN 
IV  for  the  IBM  360/67  of  the  U.   S.    Naval  Postgraduate  School.     A  copy 
of  the  program  is  contained  in  Appendix  I. 

Before  beginning  the  model  description,    it  is  desirable  to  define 
a  few  terms  commonly  used  in  the  model. 

A  replication  is  defined  as  one  complete  ASW  action.      The  action 
begins  with  the  initial  courses,    speeds,    and  locations  set  by  the  user 
for  the  ships  and  submarine  and  proceeds  until  termination  of  the 
action. 

A  run  is  a  collection  of  one  or  more  replications,    each  of  which 
has  the  same  initial  values  of  initial  courses,    speeds,    and  locations 
of  ships  and  submarine.      The  variation  among  replications  in  any  run 
are  caused  by  the  randomly  generated  occurances  within  each  replica- 
tion. 

A  high  value  group  is  a  group  of  one  or  more  ships  to  be  protected 
by  escort  ships  from  submarine  attack. 

A  screen  is  defined  as  a  barrier  formed  by  a  group  of  escort  ships, 
distributed  about  the  high  value  group,   whose  purpose  is  to  protect  the 
high  value  group  from  submarine  attack.     In  both  actuality  and  the 
model  there  are  two  general  screen  types,    the  random  patrol  area 


screen  and  the  fixed  station  screen.     In  the  random  patrol  area  screen, 
screening  ships  are  assigned  a  circular  patrol  area  with  a  specified 
radius  and  a  center  that  is  fixed  at  a  relative  range  and  bearing  from 
the  high  value  group.     The  screening  ships  randomly  patrol  this  area 
at  some  speed  greater  than  that  of  the  high  value  group.     In  the  fixed 
station  screen,    the  screening  ships  are  assigned  a  fixed  station  relative 
to  the  high  value  group  and  maintain  this   station  throughout  the 
replication. 

CPA  is  defined  as  the  closest  point  of  approach  of  a  submarine  to 
a  screening  ship. 

Detection  is  defined  as  occurring  when  a  penetrating  submarine 
reaches  a  "point  within  the  detection  range  of  a  screening  ship. 

A  formation  is  defined  as  the  group  of  ships  that  make  up  both  the 
high  value  group  and  its  screen.     The  center  of  a  formation  is  referred 
to  as   ZZ  and  may  or  may  not  be  occupied  by  a  ship. 

The  model  is  designed  to  simulate  an  ASW  environment.      The 
scenario  consists  of  a  formation,   proceeding  along  a  predetermined 
track  and  being  attacked  by  submarines  attempting  to  penetrate  the 
screen  in  specific  places.     The  entire  ASW  action  takes  place  suffi- 
ciently far  in  the  first  quadrant  of  an  X-Y  grid  to  ensure  positive 
locations  at  all  times  in  the  model.     The  ability  of  the  individual 
screen  ship  to  detect  the  penetrating  submarines  is  determined  by  a 
detection  range  that  is  uniquely  determined  for  each  screen  ship  in 
each  replication.     If  a  random  patrol  area  screen  is  used,    the  ability 


of  a  screen  ship  to  detect  a  submarine  is  also  influenced  by  the  location 
of  the  screen  ships  within  their  patrol  area;  i..  e.. ,.  at  the  extremities  of 
the  area,    the  screening  ship  may  be  several  thousand  yards  closer  or 
farther  aweiy  from  a  penetrating  submarine  than,  is  the  center  of  the 
ship's  patrol  area.     To  use  the  model,    the  user  must  supply  certain 
data  consisting  primarily  of  the  initial  locations,    courses  and  speeds 
of  all  ships,    and  some  curve  of  detection  range  v&..  probability  of 
detection  for  the  individual  ship's   sensors.      The  specific  data  are 
placed  on  control  cards  in  a  standard  format  and  a.  new  set  of  control 
cards  is  used  in  each  run.      The  flow  of  simulation- log  is   contained  in 
Figure   1. 

There  are  five  events  in  the  model.      These  are: 

1.  A  course  change  by  a  screening  ship.      This  event,    called 
whenever  a  screening  ship  reaches  the  edge  of  its  patrol  area,    gen- 
erates a  new  course  for  the  screen  ship  to  continue  its  random  patrol- 
ling. 

2.  A  detection  by  a  screening  ship.     In  this  event  the  location  of 
the  detecting  screen  ship,   detected  submarine,   and  ZZ  are  noted  at 
the  time  of  detection. 

3.  A  CPA  by  a  screening  ship.     At  CPA,    the  time;  CPA  range; 
screen  ship  and  submarine  identitites;  and  locations  of  screen  ship, 
submarine,    and  ZZ  are  stored  in  a  vector  for  reproduction  at  the  end 
of  the  run. 
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Figure   1:     Simulation  Logic   Flow 


4.  An  attack  by  the  submarine  on  a  high  value  group  ship.     The 
time  and  location  of  the  submarine  and  attacked  ship  for  each  attack 
are  noted  in  this  event. 

5.  Termination  of  the  replication.     This  event  locates  all  screen- 
ing ships,    submarines,   and  high  value  group  ships  at  replication  ter- 
mination.    The  mode  of  termination,    either  by  successful  submarine 
penetration  or  by  a  sufficient  number  of  screen  ship  detections  on.  the 
submarine,    is  indicated. 

The  measures  of  effectiveness  for  a  specified  formation  computed 
in  the  model  are: 

1.  The  number  of  attacks  on  high  value  group  ships  made  by 
penetrating  submarines. 

2.  The  number  of  times  the  penetrating  submarines  are  detected 
by  screen  ships. 

3.  The  number  of  successful  penetrations  made  by  the  submarine. 
These  measures  may  be  used  to  measure  the  effectiveness  of  a 

screen.      The  expression  "effectiveness  of  the  screen"  refers  to  the 
probability  of  the  screen  ships  detecting  the  penetrating  submarines 
and,   at  the  same  time,    detecting  the  submarine  sufficiently  far  from 
ZZ  to  prevent  attacks  on  the  high  value  group  ships.      Clearly,    a  screen 
may  be  designed  that  will  maximize  the  probability  of  detection  of  the 
penetrating  submarine;  but,    with  the  necessity  of  minimizing  the 
number  of  attacks  made  on  high  value  group  ships  as  an  added  con- 
straint,   the  problem  of  designing  an  effective  screen  is  complex. 
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With  the  measures  of  effectiveness  computed  in  the  model,    the  user  may- 
evaluate  the  relative  effect  of  changes  in  parameters;  e.g.,    changing 
the  radius  of  the  random  patrol  area  on  the  effectiveness  of  the  screen. 

It  should  be  noted  that  the  model  does  not  consider  the  aspect  of 
screen  ship  attacks  on  detected  submarines.     There  are  two  reasons 
for  this.      The  first  reason  is  that  the  inclusion  of  screen  ship  attacks 
on  detected  submarines  is  not  considered  germaine  to  the  measurement 
of  screen  effectiveness  as  previously  defined.     The  second  reason  is 
that  the  inclusion  of  screen  ship  attacks  would  greatly  complicate  the 
model,    adding  to  both  time  and  size  requirements. 
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HE.     LIMITATIONS  OF  THE  MODEL 

The  limitations  of  the  model  fall  into  the  areas  of  dimension  limita. 
tions  and  ship  maneuvering  limitations.     A  more  complete  description 
of  these  limitations  is  contained  in  the  following  paragraphs. 

A.  DIMENSION  LIMITATIONS 

The  model  is  limited  to  certain  maximum  numbers  of  ships  and 
replications  in  each  run  to  preclude  using  excessive  computer  storage. 
These  limitations  are: 

1.  Screen  ships  10 

2.  Submarines  10 

3.  High  value  group  ships  10 

4.  Replications  20 

B.  MANEUVERING  LIMITATIONS 

1.  Screening  ships  are  constrained  at  all  times  to  remain  within 
their  assigned  screening  stations.     In  the  event  of  a  detection  of  a 
submarine,    the  screen  ships  do  not  attack  but  remain  in  their  stations. 

2.  Submarine  maneuvering  is  confined  to  determination  of  a 
course  to  pass  under  ZZ. 

3.  High  value  group  ships  do  not  maneuver  but  maintain  course 
and  speed  throughout  the  replication. 
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IV.     ASSUMPTIONS 

In  the  model  certain  assumptions  have  been  made  to  reduce  the 
complexity  of  the  simulation.     These  assumptions  are  discussed  in  the 
following  paragraphs  with  respect  to  screening  ships,    penetrating 
submarines,   and  high  value  group  ships.     These  assumptions  are  con- 
sistent with  the  philosophy  of  model  design  and  the  measures  of  effec- 
tiveness derived  for  the  model. 

A.      SCREENING  SHIP  ASSUMPTIONS 

1.  Detection  ranges  are  assumed  to  be  Boolien  in  nature;  detection 
occurs  with  probability  one  when  the  penetrating  submarine  is  within 
this  range,    and  detection  does  not  occur  otherwise..     The  detection 
ranges  are  generated  uniquely  for  each  ship  in  each  replication. 
Detection  ranges  are  generated  at  the  beginning  of  each  replication 
using  the  following  procedure:     The  model  user  provides  some  curve 

of  detection  probability  vs.    detection  range  for  each  sensor  that  is 
associated  with  the  screen  ships.     A  number  is  randomly  drawn  from 
a  uniform  distribution  on  [0,  1];  this  number  is  the  detection  probability 
for  the  ship  in  that  replication.      The  range  corresponding  to  this  detec- 
tion probability  for  the  sensor  is  the  detection  range  of  the  ship  for  that 
replication. 

2.  As  indicated  above,    the  model  is  designed  making  an  assumption 
of  a  uniform  distribution  on  [0,  l]  for  detection  probability.     If  the  user 
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feels  that  detection  probabilities  are  distributed  according  to  some  dis- 
tribution other  than  uniform,   the  model  subroutine  generating  screen 
ship  detection  ranges  must  be  rewritten  by  the  user  to  reflect  this 
different  distribution. 

3.  Screening  ships  remain  within  their  assigned  stations  at  all 
times  in  the  model. 

4.  Screen  ships  maintain  a  constant  speed  throughout  each  run. 

5.  Upon  detection  of  a  penetrating  submarine,,  the  detecting  ship 
does  not  attack  the  submarine.     The  detection  of  the  submarine  is  noted, 
and  the  replication  continues. 

6.  Once  a  screening  ship  has  gained  contact  on  a  penetrating 
submarine,    the  interactions  of  that  screen  ship  and  that  submarine 
are  terminated;  i.e.,   as  the  replication  continues,    the  screen  ship 
does  not  redetect  the  submarine  or  calculate  a  new  CPA  for  the  sub- 
marine with  each  new  random  course  change  of  the  screen  ship  in  its 
patrol  area. 

B.      SUBMARINE  ASSUMPTIONS 

1.       The  penetrating  submarine  conducts  no  evasive  maneuvering 
in  the  model.      This  assumption  corresponds  to  actual  submarine 
maneuvering  practices.      The  submarine  fire  control  solution  for  a 
penetration  yields  a  course  to  pass  under  ZZ..     Elaborate  submarine 
maneuvering  in  a  screen  penetration  is  avoided,    if  possible,    as  any 
maneuvering  has  a  detrimental  effect  on  the  submarine  fire  control 
solution. 
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2.       The  penetrating  submarine  conducts  no  evasive  maneuvering 
while  penetrating  the  screen,    nor  does  it  attack  screen  ships.     The 
submarine  does  conduct  attacks  on  high  value  group  ships  when  arriving 
within  a  weapon's   release  range  set  by  the  user  for  each  run.      These 
attacks  allow  records  of  the  effects  of  submarine  penetrations  to  be 
maintained. 

C.      HIGH  VALUE  GROUP  ASSUMPTIONS 

1.  The  high  value  group  ships  do  not  maneuver  in  the  model  but 
maintain  the  initial  course  and  speeds  assigned  by  the  user  for  each 
run. 

2.  High  value  group  ship  number   1  is  always  located  at  ZZ.     If 
the  formation  for  the  high  value  group  does-  not  have  a  ship  at  ZZ,    high 
value  group  ship  number  1  is  a  dummy  at  ZZ,    and  the  maximum  high 
value  group  size  is  reduced  to  nine. 
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V.     INPUTS 

Inputs  to  the  model  are  by  control  cards,    one  set  for  each  run. 
The  entries  on  the  control  cards  are  all  right  adjusted;  i.  e.  ,    the  entry- 
values  set  as  far  to  the  right  as  possible  in  the  column  width  specified. 
The  input  cards  are  ordered  sequentially- as:  listed  below  and  follow  the 
source  deck  directly. 

A.      CARD  NUMBER  1 

The  entry  values  of  control  card  1  are  integer  values. 
Columns  Definition  Range 


1-4  Number  of  screen  ships  in  run  l-N-10 

5-8  Number  of  submarines  in  run  l-N-10 

9-12  Number  of  high  value  group  ships  in  run  l-N-10 

13-16  Number  of  replications  in  run  l-N-20 

17-20  Number  of  detections  of  submarine  by  screen 

ships  to  terminate  replication  1-N-ll 

21-24  Print  screen  ship  detection  range  for  each 

replication  &c  average  detection  range  for 

each  ship  for  the  run  0,No;l,Yes 

25-28  Print  out  CPA  &  locations  of  each  ship  & 

submarine  in  replication  0,No;l,  Yes 

29-32  Print  out  locations  of  screen  ship,     sub- 

marine,   &  ZZ  at  time  of  a  penetrating  sub- 
marine detection  by  a  screen  ship  0,  No;  1,  Yes 
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Columns  Definition  Range 

33-36  Minimum  submarine  penetration  speed  0-N 

37-40  Print  out  locations  of  screen  ships,    sub- 

marines,   &  ZZ  at  termination  of  each. 

replication  0,No;l,Yes 

41-44  Print  out  event  calendar  for  each  replication       0,No;l,  Yes 

45-48  Print  out  screening  ship  information  for  each 

replication  0,No;l,Yes 

49-52  Print  out  high  value  group  ship  information 

for  each  replication  0,No;l,  Yes 

53-56  Print  out  submarine  information  for  each 

replication  0,No;l,Yes 

57-60  Submarine  weapon  release  range  (100's  of 

yards) 
61-64  Print  out  locations  of  high  value  group  ship  & 

submarine  at  weapon  release  for  each 

replication  0,No;l,Yes; 

B.      CARDS  FOR  SCREEN  SHIPS 

Cards  for  screen  ships  follow  sequentially  the  first  control  card. 
These  entries  and  those  of  following  cards  must  be  non-negative  and 
may  be  integer  or  floating  point.      For  each  screen  ship,    there  is  one 
card  with  the  following  format: 
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Columns  Definition 


1-8  Screen  ship  number 

9-16  Initial  course  (in  radians) 

17-24  Ship  speed 

25-32  Initial  X-position 

33-40  Initial  Y-position 

41-48  Sensor  type 

57-64  Patrol  station  radius 

65-72  Initial  X_position  of  center  of  random  patrol  area 

73-80  Initial  Y-position  of  center  of  random  patrol  area 

C.  CARDS  FOR  PENETRATING  SUBMARINES 

Cards  for  penetrating  submarines  follow  directly  after  the  cards 
for  the  screen  ships.      For  each  submarine,    there  is  one  card  with 
the  following  format: 

Columns  Definition 

1-8  Submarine  number 

17-24  Speed  (    -    minimum  penetration  speed) 

25-32  Initial  X-position 

33-40  Initial  Y-position 

D.  CARDS  FOR  HIGH  VALUE  GROUP  SHIPS 

Cards  for  the  high  value  group  ships  follow  directly  after  the 
cards  for  the  penetrating  submarines.      For  each  high  value  group 
ship,    there  is  one  card  with  the  following  format: 
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Columns  Definition 


1-8  High  value  group  ship  number 

9    -16  Course 

17-24  Speed 

25-32  Initial  X-position 

33-40  Initial  Y-position 
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VI.     OUTPUTS 

The  outputs  of  the  model  are  of  two  types,    general  statistical  out- 
puts and  specific  format  outputs.     The  general  statistical  outputs  are 
always  provided  by  the  model,   while  the  specific  format  outputs  are 
provided  to  the  user  only  if  they  are  desired^.     To  have  specific  format 
outputs,    the  user  codes  control  card  1  in  the  appropriate  column  with 
a  1  for  each  specific  format  output  desired  and  a  0  for  those  specific 
format  outputs  not  desired.     Examples  of  general  statistical  and 
specific  format  outputs  are  included  in  Figures  2,    3,    and  4. 

A.      GENERAL' STATISTICAL  OUTPUTS 

1.  The  number  of  attacks  on  high  value  group  ships  by  sub- 
marines for  each  replication  and  the  total  for  the  run. 

2.  The  number  of  detections  made  on  the  penetrating  submarine 
by  screening  ships  during  each  replication.      The  upper  limit  of  this 
value  is  set  on  control  card  one  for  each  run. 

3.  The  number  of  successful  penetrations  to  ZZ  made  by  penetra- 
ting submarines  in  each  replication  and  the  total  for  the  run. 

4.  Average  detection  range  of  the  penetrating  submarine  from  ZZ 
for  each  replication;  average  detection  range  and  its  standard  deviation 
of  the  penetrating  submarine  from  ZZ  for  each  run. 
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B.      SPECIFIC  FORMAT  OUTPUTS 

1.  Individual  screen  shiprs  generated  detection  ranges  for  each 
replication  and  the  average  generated  detection  range  for  the  run. 

2.  CPA  of  each  screen  ship  to  a  submarine  in  each  replication. 

3.  Locations  of  ZZ,    screen  ship  making  detection,    and  the 
detected  submarine  at  time  of  detection  in  each  replication. 

4.  Locations  of  ZZ,    screen  ships,    and  submarines  at  the  ter- 
mination of  each  replication. 

5.  Event  calendar  for  each  replication.. 

6.  Screen  ship  information  at  termination  of  each  replication. 
This  information  consists  of: 

a.  Screen  ship  number 

b.  Course  (in  radians) 

c.  Speed 

d.  Ship  X-position 

e.  Ship  Y-position 

f.  Sensor  type 

g.  Generated  detection  range 
h.       Patrol  station  radius 

i.        X-position  of  center  of  random  patrol  area 

j.        Y-position  of  center  of  random  patrol  area 

k.       Time  of  next  random  course  change 

1.        Detection  status  of  ship  (1  if  ship  made  detection, 

0  otherwise) 
m.      Time  of  last  position  update 
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7.  High  value  group  ship  information  at  termination  of  each 
replication.      This  information  consists  of: 

a.  High  value  group  ship  number 

b.  Course  (in  radians) 

c.  Speed 

d.  X-position 

e.  Y-position 

f.  Time  of  last  position  update 

g.  Time  first  submarine  passes  under  ZZ 

8.  Submarine  information  at  termination  of  each  replication. 
This  information  consists  of: 

a.  Submarine  number 

b.  Course  (in  radians) 

c.  Speed 

d.  X_position 

e.  Y_position 

f.  Time  of  last  position  update 

g.  Number  of  times  the  submarine  has  been  detected  by- 
different  screen  ships.     The  upper  limit  of  this  value  is  set  on  control 
card  1  for  each  run. 

9.  Locations  of  submarine  and  target  ship  at  weapon  release  in 
each  replication. 
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1.       Examples  of  General  Statistical  Outputs 
a.       Output  for  each  replication 

WEPLAUNCH  DETECTIONS  PENETRATION 


b.       Output  for  each  run 

PENETRATIONS         WEPLAUNCH         AVTJET        STDDEV 

8  12  6057.6  2386.65 

AVERAGE  RANGE  OF  SUBMARINE  DETECTION  FROM  ZZ 
IN  EACH  REPLICATION 

7765.18         8809.04         0.0         4286.71         8633..73         6465.39 

2.       Examples  of  Specific  Format  Outputs 

a.  Individual  screen  ship's  generated  detection  ranges  for  each 
replication  &t  the  average  generated  detection  range  for  the  run 
SHIP      AVRANGE      REP  1      REP  2      REP  3      REP  4      REP  5 

1  3049.9  0.0         1097.8       0.0  3380.6     7427.6 

b.  CPA  of  each  screen  ship  to  a  submarine  in  a  replication 

SHIP     SUB      TIME      CPA      SHIPXLOC      SHIPYLOC      SUB  XLOC 

3.0        1.0  73.        7245.  44238.  68293.  37006. 

SUB  YLOC      ZZ  X  LOC      ZZ  Y  LOC 
67857.  40000.  64319. 

c.  Location  of  ZZ,    screen  ship  &  detected  submarine  at  time  of 

detection 

DETECTION  TIME  64  BY  DD  4  ON  SUBMARINE  1  AT 
X=36733.00    Y-62409.00 

Figure  2:     Examples  of  Outputs  Available  in  the  Model 
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d.       Locations  of  ZZ,    screen  ships,    &  submarines  at  replication 
termination 


TIME 


86 


TIME 


SCREEN  SHIP 


86 


TIME 


MAIN  BODY  SHIP 


SUBMARINE 


86 


X-POSIT 

42180.0 

X-POSIT 

40000.0 

X-POSIT 

39910.0 


Y-POSIT 

72546.0 

Y-POSIT 

68652.0 

Y-POSIT 

68724.0 


e.       Event  calendar  for  each  replication 


TIME 


EVENT  NO. 


SHIP  SUB 


MAIN  BODY 


86 


f.        Screen  ship  information  at  termination  of  each  replication 
SHIP      COURSE      SPD      X- COORD      Y-COORD      SONAR 


1.62 


15. 


35434. 


72029. 


0. 


DETRANGE       PATRAD      PATX        PATY      CCTIME      DS 


4691.3 


1800.    35830.   71074. 


90. 


XYTIME 


86. 
g.       High  value  group  ship  information  at  the  termination  of  each 
replication 
SHIP      COURSE      SPD      X_  COORD      Y-COORD      XYTIME 


1 


0.0 


10.  40000. 


68652. 


86.0 


EVT9TIME 

86. 
Figure  3:     Examples  of  Outputs  Available  in  the  Model 
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h.       Submarine  information  at  the  termination  of  each  replication 
SHIP        COURSE        SPD        X- COORD        Y- COORD        XYTIME 

1.  1.28  7.0  39910.  68724.  86.0 

TIMES  DETECTED 
1. 
i.        Locations  of  submarine  h  target  ship  at  -weapons'   release  in 
each  replication 

TIME        SUB  NO        SUB  X        SUB  Y        TARGET  NO        TAR  X 
74  1  37230.        67927.  1  40000. 

TAR  Y         LAUNCH  RANGE 
64654.  4000. 


Figure  4:     Examples  of  Outputs  Available  in  the  Model 
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VII.     MODEL  USE 

As  an  example  of  the  use  of  the  ASW  simulation  model,   the  prob- 
lem of  detection  for  a  screen  of  four  ships  is  optimized  against  a 
penetrating  submarine  threat.      The  use  of  the  term  optimization  in 
the  context  of  this   section  refers  to  a  search  for  an. improved,    though 
not  necessarily  optimal,    solution.     The  screen  is  specified  to  be  a 
random  patrol  area  screen  -with  the  four  screen  ships  symmetrically 
distributed  about  the  formation  base  course  in  patroL  areas  whose 
centers  are  equally  spaced.      This  symmetric  distribution  of  screen 
ships  will  cover  a  specific  arc  in  front  of  the  high  value  group.      The 
screen  ships  must  be  located  a  distance  from  ZZ  that  is  at  least 
equal  to  maximum  submarine  -weapons'  launch  range.      The  submarine 
threat  is  known,    with  the  submarine  specified  to  pass  between  screen 
stations.      The  size  of  the  main  body  is  unspecified.      (See  Figure  5) 

A.      PROCEDURE 

The  procedure  used  for  optimization  in  this  thesis  is  that  of 
response  surface  methodology  (RSM).     This  procedure  focuses  on 
the  class  of  simulations  for  which  a  single  output  variable  exists, 
and  the  input  and  output  variables  are  quantitative  and  approximately 
continuous. 
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Figure  5:     Four  Ship  Random  Patrol  Area  Screen  with 
Specified  Submarine  Penetration  Tracks 
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If  there  exists  a  functional  relationship  between  the  quantitative 
output  variable  and  the  quantitative  input  variables,    this  relationship 
may  be  expressed  as: 

y  =  g  (xj,.. .,xk) 
■where:  y    denotes  the  quantitative  response  variable 

x.  denotes  the  it£L  input  variable 
g    denotes  the  functional  relationship 
and  the  value  of  the  response  variable  y  lies  on  some  k  +  1  dimen- 
sional response  surface. 

Assuming  this  functional  relationship  g  to  be  continuous  and 
differentiate,    we  will  assume  that  a  Ta.ylor  series  expansion  of  g 

can  be  expressed  as: 

k  K       k 

y  -  g(X)   =  a     +    Z    aiXi    +  2L  II  aijXiXj    +  •  •  • 

In  the  above  notation,     y   is  a  measure  of  screen  effectiveness  deter- 
mined for  a  screen  by  an  average  of  the  number  of  submarine  detec- 
tions over  all  replications  in  the  runs;    and  X.  is  the  value  of  the  i^*1 
input  parameter  for  the  screen. 

The  problem  using  the  simulation  modeL  is  that  we  do  not  know 
the  functional  relationship  g,  yet  we  wish  to  maximize  g(X)  subject 
to  some  set  constraints;    i.  e.  ,    maximize  y  subject  to  maximum  and 

minimum  values  for  X..      The  model  is  used  to  estimate  the  coefficients 

1 

of  the  above  Taylor  series  expansion  by  sampling.      That  is,    the  model 
will  be  used  to  obtain  values  for  y  for  specified  values  of  X.;  and  the 
technique  of  least  squares  can  be  used  to  estimate  the  values  of  the 
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coefficients  of  the  Taylor  series  expansion.     Having  an  approximation 

for  the  functional  form  of  g(X),     it  is  possible,    from  this  functional 

form,    to  estimate  the  gradient  of  g  and  to  use  this  estimated  gradient 

to  determine  the  direction  to  move  for  further  sampling  in  an  effort 

to  maximize  g(X). 

If  the  initial  samples  used  to  determine  the  values  of  the  estimators 

of  the  coefficients  of  the  Taylor  series  expansion  are  taken  from  a 

rather  restricted  area,    it  is   reasonable  to  make  the  assumption  that 

g  is  linear  throughout  this  area.      That  is,    the  functional  form  of  g(X) 

may  be  written: 

k 

y  =  g(X)  =  a0  +   £   aixi 

and  that  the  gradient  of  g  is: 

(a -.,...  ,  a    ). 

The  sampling  procedure  is  to  observe  at  least  k  +  1  values  of  y  at 

k  +  1  different  values  of  X  and  to  estimate  the  values  of  (a,,  a i .....  a.  ) 

—  U       -1  k 

by  least  squares.     A  large  number  of  simulation  runs  with  different 

values  of  the  X-'s  may  be  made  to  determine  the  values  of  the  a.'s  by 

least  squares.      By  considering  only  the  largest  and  smallest  values  of 

the  X-'s  within  the  restricted  area,    there  are  2      experiments  that  may 

be  performed  without  redundancy  to  determine  the  value  of  the  a.'s. 

l 

The  minimum  number  of  experiments  that  will  uniquely  determine  the 
values  of  the  a-'s  in  k  +  1  dimensional  space  is  k  +  1.  Due  to  the  high 
cost  of  computer  time,    the  author  did  not  consider  it  necessary  to  use 
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more  than  k  +  different  values  of  the  input  variables  to  determine 

the  value  of  the  a.'s. 

1 

To  ease  the  computational  complexity  of  the  solution,    the  user  may- 
code  each  X.  with  +1  for  its  largest  value.1  and  -L  for  its  smallest  value. 
By  coding  the  input  variables  as   +1,    the  assumed  linear  function 
relationship  g  is  transformed  to  a  new  linear  functional  relationship 

of  the  form: 

k 

y         =  b     +    Z  b.x. 
0  i       i   i 

where        x.        =  X^     _     A  + 


A 


=    maximum  X^  +  minimum  X^ 


2 


/\     =  maximum  X.    _  minimum  X. 

i  l 


and  bQ      =  aQ  +  £  a.    £,+ 

b.       =  a.    A" 

By  using  k  +  1  experiments,    a  unique  solution  for  the  estimated  values 

of  b     through  b,    is  obtained  by  the  simultaneous   solution  of  k  +  1 
Ok  ' 

linear  equations. 

The  estimates  of  b     through  b     obtained  by  sampling  are  then 
used  as   estimates  of  the  gradient  of  the  response  surface.      By  using 

the  gradient  as  a  guide,    one  expects  to  proceed  to  an  improved  solution 

k 
V    =    b     +  w  b.  2 

°  1        l 

where  y  denotes  the  predicted  response  along  the  path 
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(wb1>  wb2,  .  .  .  ,wbk) 
where  wis  a  positive  constant  to  indicate  we  are  moving  in  the  direction 
of  increasing  values  of  y.     That  is,   we  shall  observe  a  new  value  of  y 
at  the  point  (wb    ,  wb^,  .  »  .  wb    )  and  hope  that  the  observed  value  of  y 
will  be  close  to  the  predicted  value  of  y  using  y  above.      It  should  be 
noted  that  with  increasing  w,    the  accuracy  of  the  predicted  response 
may  become  unreliable  unless  the  unknown  functional  relationship  is, 
in  fact,   linear;  and,    thus,    this  path  should  be  followed  with  caution. 
As  long  as  observed  values  for  y  continue  to  increase  along  the  path, 
one  can  continue.     If  the  linear  approximation  for  the  unknown  func- 
tional relationship  has  seemed  satisfactory  along  the  path,    the  appear- 
ance of  stationary  or  decreasing  values  for  y  as  the  path  is  continued 
to  be  followed  indicates  that  the  linear  approximation  and  gradient  is 
calculated.      The  new  path  indicated  by  this  gradient  is  followed  until 
either  the  optimal  solution  is  reached,    or  a  constraint  is  reached;  or 
a  new  linear  approximation  and  gradient  must  be  calculated. 

B.      THE  PROBLEM 

For  the  problem  being  considered,    the  screen  is  specified  to  be  a 
random  patrol  area  screen  with  the  ship  located  on  bearing  300    ,    340    , 
020    ,    and  060      relative  to  ZZ  respectively.      The  screen  center  of 
each  patrol  area  is  to  be  no  greater  than  10,000  yards  from  ZZ  and 
no  less  than  4,  000  yards  from  ZZ.     The  submarine  is  to  penetrate 
from  320    ,    000    ,    and  040      relative  to  ZZ.      The  number  of  ships  in  the 
high  value  group  is  unspecified.      (See  Figure  5) 

31 


The  model  provides  for  each  run  the  number  of  successful  pene- 
trations made  by  the  submarine.     Therefore,   using  the  number  of 
unsuccessful  penetrations  as  a  measure  of  screen  effectiveness  to 
optimize  using  RSM,    one  must  first  consider  what  variables  in  the 
model  affect  this  measure  of  effectiveness.     The  variables  that  affect 
the  output  parameter  of  effectiveness  y  are:     ship  speed,,  radius  of 
random  patrol  area,    the  range  of  the  center  of  the  random  patrol  area 
from  ZZ,    ship  detection  range  and  the  location  of  the  ship  within  its 
random  patrol  area.     The  first  three  variables,   ship  speed,    radius  of 
patrol  area,    and  the  range  of  the  center  of  the  patrol  area  from  ZZ, 
are  variables  over  which  the  user  has  direct  control  in  the  model. 
Ship  detection  range  and  the  location  of  the  ship  within  its  random 
patrol  area  are  random  within  the  model  and,    therefore,    not  directly 
controlled. 

For  this  problem,    it  was  decided  to  limit  the  number  of  computer 
runs  to  eleven,  with  ten  replications  in  each  run,  and  to  optimize  the 
screen  effectiveness  using  no  more  than  those  runs.     Screen  ship 
speed  was  fixed  at  15  knots  with  the  speed  of  the  high  value  group  at 
10  knots.      Thus,    there  are  two  variables  to  consider  in  optimizing  the 
screen:     radius  of  the  patrol  area  (X,)  and  the  range  of  the  center  of 
the  patrol  area  from  ZZ  (X2). 

For  sampling,  the  initial  area  considered  was  bounded  by  the 
values  of  1,  000  yards  and  3,  000  yards  for  X,  and  8,000  yards  and 
10,  000  yards  for  X^.      The  coded  variables  are: 
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X 
x     =  1     -     2,000 


1  1,000 

x     _         X2    -    9,000 


2  1,000 

The  initial  sample  points  and  their  associated  responses  were: 
Sample  point  X  X_  x.  x->  y 

1  1,000  10,000  -1  I  0L.533 

2  1,000  8,000  -I  _I  0.567 

3  3,000  8,000  -1  _L  0.667 
The  solution  of  the  three  simultaneous  linear  equations  yielded 

the  following  linear  approximation: 

y  =  0.457  +  0.017x     -  0. 075x2 
Since  the  largest  observed  value  of  the  response  was  0.667,    it 
was  decided  to  choose  w  such  that  the  predicted  response  would  be  in 
the  neighborhood  of  0.700.      The  w  used  was  41,    which  corresponded 
to  y  =  0.697.     We,   therefore,    sampled  y  at  the  point 

(41  (0.017)  ,    41  (-0.075)) 
The  data  for  the  next  observation  was  as  follows: 
Sample  point  w  X,  X  x  x  y  y 

4  41        2,700      5925     w(0.017)    w(-0.075)    0.667    0.697 
The  value  of  the  observed  response  for  this  run  was  0. 667,   a  value 
less  than  predicted  but  still  not  less  than  the  greatest  previously 
observed  value.     It  was  decided  to  continue  along  this  path  to  see  if 
the  solution  would  improve.     Accordingly,    a  value  for  w  of  49, 
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corresponding  to  a  predicted  response  of  0.747,   was  chosen.     We 
sampled  y  at  the  point 

(49  (0.017)  ,   49  (-0.075)).. 
The  data  for    this  next  observation  was  ata  fallows: 
Sample  point        w         X  X?  x  x-  y  y 

5  49        2800       5325      49(Q..017)      49(-0.075)      0.767      0.747 

By  examination  of  the  differences  between  predicted  and  response  out- 
put for  the  previous   runs,    it  appeared  that  the  linear  approximation 
was  not  valid.     Nevertheless,    as  long  as  the  response  continued  to 
increase,    the  path  indicated  by  our  approximation  of  the  gradient  of  g 
appeared  to  indicate  the  direction  to  an  improved  solution.      Continuing 
along  the  path,    we  next  selected  w  =  58,    a  value  corresponding  to 
y  =  0.  800,    and  sampled  at  the  point 

(58  (0.017)  ,    58  (-0.017)). 
The  data  for  this  next  observation  was  as  follows: 

A 

Sample  point        w       X  X_  x.  x  y  y 

6  58       3000       4650       58(0.017)       58(-0.075)        0.667      0.800 

As  the  value  of  the  response  has  decreased,    it  was  necessary  to 
return  to  the  point  of  maximum  response  and  calculate  a  new  linear 
relationship  and  gradient.     Again,    three  sample  points  were  used  with 
the  previous  highest  response  point  as  one  of  these.      These  sample 
points  and  their  associated  responses  were: 
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Sample  point  X.  X  x  x  v 

1  2  L  2 

7  1800  5325  -I  1  0.733 

8  2800  5325  I  1  0.733 

9  2800  4325  I  -1  0.800 

The  solution  of  the  simultaneous  linear  equations  associated  with 
these  sample  points  yielded  the  following  linear  approximation: 

y  =  0.7665  +  Ox,    -    0. 0335x   .. 

J  1  2 

It  is  of  interest  to  note  that  the  value  of  the  estimator  correspond- 
ing to  the  value  of  the  radius  of  the  patroLarea  is  zero.      This  indicates 
that  the  value  of  the  radius  is  at  its  optimal  value  and  can,    therefore, 
no  longer  contribute  to  an  improved  solution.     A  value  of  x     =0  cor- 
responding to  X,    =  2,300  yards  was,    therefore,    assigned  in  the  final 
runs  of  the  model. 

We  decided  to  choose  w  =  30  such  that  the  predicted  response 
would  be  y  -  0.  800.     We  sampled  y  at  the  point 

(30    (0),    30  (-0.0335)). 
The  data  for  this  next  observation  was  as  follows: 
Sample  point  X  X  x,  x  y  y 

10  2300        4820  0        30(-0.0335)        0.800        0.800 

It  was  realized  that  the  predicted  response  of  sample  point  10  was  not 
an  improvement  over  the  best  previous  value;   but  the  value  of  X_  was 
approaching  the  overall  minimum  constraint  of  4,  000  yards,    and  a 
more  cautious  step  was  taken. 
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Continuing  along  the  path  for  the  final  run,    a  value  of  w  =  74 
corresponding  to  y  =  0.  850  was  chosen.     We  s-ampled  y  at  the  point 

(74  (0),    74  (-0.0335)). 
The  data  for  this  final  observation  was  as  foLLows: 
Sample  point  X  X  x  x  y  y 

11  2300  4085  0  74(-Q..G335)        0.833        0.850 

With  the  completion  of  the  eleven  runs  of  the  model,    the  optimal 
screen  for  the  given  constraints  was  found  to  have  values  of  X,   and 
X_  of  2,300  yards  and  4,085  yards   respectively..     The  optimal  response 
for  these  values  of  X     and  X     was  y  =  0.  833;  this  means  that  the  op- 
timized  screen  can  be  expected  to  detect  a  penetrating  submarine  prior 
to  the  submarine  attacking  the  high  value  group  83.  3%  of  the  time. 
With  additional  runs,    the  response  could  possibly  be  improved;  but, 
with  only  the  variable  X?  available  to  vary  and  with  its  value  so  close 
to  its  lower  limit  of  4,000  yards,    the  improvement  would  probably  be 
marginal.     A  summary  of  the  above  steps  is  contained  in  Figure  6. 

In  this  sample  example,    a  simple  screen  that  had  several  con- 
straints was  considered.      By  using  RSM,    the  screen  variables  were 
optimized  against  the  threat  of  submarine  penetration.      Using  the 
same  techniques,    the  model  may  be  used  to  optimize  more  complex 
screens  with  additional  variables. 
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Run  No. 

Xl 

X2 

Xl 

X2 

A 

y 

1 

1,000 

10,000 

.1 

1 

2 

1,000 

8,000 

_1 

-1 

-- 

3 

3,000 

8,  000 

1 

-1 

Y         Comments 

Q.533     Value  for  path  of 

steepest  ascent  are: 

0.567     (Z000  +  17w,  9000-75w) 
in  Xi ,  Xi   space 

0.667     y=0.457+w(0. 005914) 


7 
8 


2,700  5,925 

2,800  5,325 

3,000  4,650 

1,800  5,325         -1 

2,800  5,325  1 


0.697     0.667    w  =  41 

0.747     0.767    w  =  49  point  of 

maximum  response 

0.80       G.667    w  =  58 

0.733 

0.733     y=0.7665+w(0. 00112225) 
b  =0.0  in  this   solution 


2,  800 


4,  325 


1    -1    --    0.80 


10 


11 


2,300     4,820 
2,300     4,085 


0.8  0.8  w  =  30 


0.85        0.833     w  =  74 


Figure  6:    Summary  of  the  Search  for  Improved  Solution 
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VIII.      CONCLUSIONS  AND  RECOMMENDATIONS 

The  model  presented  in  this  thesis  is  a  simulation,  model  of  ASW 
interactions  between  a  formation  and  penetrating  submarines.     A  pro- 
cedure has  been  demonstrated  in  which  the  model  was  used  to  improve 
the  effectiveness  of  a  proposed  screening  tactic  against  a  known  sub- 
marine threat. 

An  attempt  to  use  the  model  to  measure  the  effectiveness  of 
various  screen  tactics  was  examined.     The  procedure  used  in  this 
attempt  was  to  determine  the  average  detection  range  of  penetrating 
submarines  from  ZZ  for  each  screen  and  to  compare  these  average 
ranges  by  the  analysis  of  variance  techniques.     This  procedure  to  rank 
screen  effectiveness  in  order  of  average  detection  range  from  ZZ, 
with  the  screen  having  the  greatest  average  detection  range  being  the 
most  effective,    had  intuitive  appeal.     Unfortunately,    the  standard 
deviations  of  the  average  detection  ranges  were  so  great,    on  the  order 
of  two  to  three  times  the  differences  between  average  detection  ranges 
of  the  screens,    that  no  statistical  interferences  could  be  made  with 
this  procedure. 

In  an  effort  to  use  the  same  general  procedure,    an  attempt  was 
made  to  compare  the  effectiveness  of  proposed  screen  tactics  by  com- 
paring the  average  time  until  first  detection  of  a  penetrating  submarine. 
This  procedure  also  yielded  no  statistically  significant  results  for 
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the  same  reason  as  the  previous  proposal.      The  results  of  these 
attempts  are  not  included  in  this  thesis. 

The  model  has  not  been  successfully  used  to  directly  compare 
the  effectiveness  of  various  proposed  screening,  tactics;  but  it  may 
be  used  as  shown  in  the  preceding  chapter  to  maximize  the  effective- 
ness of  a  screen,   with  certain  constraints,   againsta  specific  sub- 
marine threat. 
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APPENDIX  A.     COMPUTER  PROGRAM 
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APPENDIX  B:      LIST  OF  PROGRAM  VARIABILES 

This  appendix  contains  the  name  and  definition  of  the  symbols 
appearing  in  the  computer  program.  This  list  contains  only  those 
variables  using  storage  space  in  the  program  and  does  not  contain 
those  used  only  for  temporary  storage. 
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Name 


Definition 


CPA(I,  J,  K)  CPA  information  for  the  Ith  ship:  air.  the  Ith  submari 


ne 


CPA(I,  J,  1) 
CPA(I,  J,  2) 


in  a  replication.     The  values  in.  this  vector  are: 
Screen  ship  number 
Penetrating  submarine  number 


CPA(I,  J,  3)  Time  of  CPA 


CPA(I,  J,  4)  CPA  range 


CPA(I,  J,  5) 
CPA(I,  J,  6) 


Screen  ship  X  location  at  CPA 
Screen  ship  Y  location  at  CPA 


CPA(I,  J,  7)  Penetrating  submarine  X  location  at  CPA 


CPA(I,  J,  8) 
CPA(I,  J,  9) 


Penetrating  submarine  Y  location  at  CPA 


X  location  of  ZZ  at  CPA 


CPA(I,  J,  10)         Y  location  of  ZZ  at  CPA 
DECAL  (I,  J)  Dummy  event  calendar,  a  copy  of  all  events   entered  in 

the  event  calendar.      The  I       entry  in  DECAL  contains: 
DECAL(I,  1)         Event  time 
DECAL(I,  2)         Event  number 
DECAL(I,  3)         Number  of  screen  ship  in  event 


DECAL(I,  4) 
DECAL(I,  5) 
DESKIL 


Number  of  penetrating  submarine  in  event 
Number  of  main  body  ship  in  event 

An  index  of  successful  attacks  by  a  penetrating  sub. 
marine  on  a  main  body  ship  in  a  replication 
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Name 


Definition 


DSSV(I,  J) 
DSSV(I,  1) 


Dummy  submarine  vector  containing  initiaH information 
Number  of  I      penetrating  submarine 


DSSV(I,  3)  Speed  of  submarine 


DSSV(I,  4) 


DSSV(I,  5) 


Initial  X  location  of  submarine 


Initial  Y  location  of  submarine 


rth 


DSTAT(I,  J,  K)         Detection  range  from  ZZ  in  the  I       replication  for  the 

R**1  submarine  by  the  J*     screen  ship 


DSV(I,  J) 
DSV(I,  1) 
DSV(I,2) 
DSV(I,  3) 
DSV(I,  4) 
DSV(I,  5) 
DSV(I,  6) 
DSV(I,  8) 

DSV(I,  9) 
DSV(I,  10) 
DTAR(I,  J) 


DTAR(I,  1) 

DTAR(I,  2) 

DTAR(I,  3) 

DTAR(I,  4) 
DTAR(I,  5) 


Dummy  screen  ship  vector  containing  initial  information 

Number  of  I       screen  ship 

Initial  screen  ship  course  (in  radians) 

Screen  ship  speed 

Initial  screen  ship  X  location 


Initial  screen  ship  Y  location 


Screen  ship  sonar  type 

Screen  ship  patrol  area  radius 

Initial  X  location  of  center  of  assigned  patrol  area 

Initial  Y  location  of  center  of  assigned  patrol  area 

Dummy  main  body  ship  vector  containing  initial  in. 


formation 


rth 


Number  of  I      main  body  ship 

Initial  main  body  ship  course  (in  radians) 

Main  body  ship  speed 

Initial  X  location  of  main  body  ship 
Initial  Y  location  of  main  body  ship 


75 


Name 

EA1 

EA2 

EA3 

ECAL(I,  J) 
ECAL(I,  1) 
ECAL(I,  2} 
ECAL(I,  3} 
ECAL(I,  4) 
ECAL(I,  5) 

EN 

ET 

KT 

MAXEV 

MBKIL 


MBKILL 
MSTAT(I) 
MSTAT(l) 
MSTAT(2) 
MSTAT(3) 
MSTAT(4) 
MSTAT(5) 


Definition 


Number  of  a  screen  snip 

Number  of  a  penetrating  submarine 

Number  of  a  main  body  ship 

th 
I       entry  in  event  calendar 


Event  time 


Event  number 


Number  of  screen  ship  in  event 

Number  of  penetrating  submarine  in  event 


Number  of  main  body  ship  in  event 


Event  number 


Event  time 


Index  of  total  number  of  events  entered  in  DECAL(I,  J) 

Maximum  number  of  events  allowed  in  ECAL(I,  J) 

Number  of  successful  submarine  penetrations  in 

replication 

Number  of  successful  submarine  penetrations  in  run 


Control  vector  for  each  run 


Number  of  screen  ships  in  run 


Number  of  submarines  in  run 


Number  of  main  body  ships  in  run 

Number  of  replications  in  run 

Number  of  detections  of  submarine  by  screen  ships  to 


terminate  run 
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Name 


Definition 


MSTAT(6)  Print  screen  ship  detection  range  for:- each  replication 

&  average  detection  range  far  each  s.cr£en  ship  during 

run 
MSTAT(7)  Print  out  CPA  &  locations  of  each  ship  h  submarine  in 

replication 
MSTAT(8)  Print  out  locations  of  screen  ship,    submarine,    and 

ZZ  at  time  of  detection  of  a  penetrating  submarine  by 


MSTAT(9) 
MSTAT(IO) 


MSTAT(ll) 
MSTAT(12) 
MSTAT(13) 

MSTAT(14) 
MSTAT(15) 
MSTAT(16) 


MTARK 

PI 

R 


a  screen  ship 

Minimum  submarine  penetration  speed 

Print  out  locations  of  screen  ships,    submarines,    & 

ZZ  at  termination  of  each  replication 

Print  out  event  calendar  for  each  replication 

Print  out  ship  vectors  for  each  replication 

Print  out  main  body  vectors  for  each  replication 

Print  out  submarine  vectors  for  each  replication 

Submarine  weapon  release  range  (100's  of  yards) 

Print  out  locations  of  main  body  ship  &c  submarine  at 

weapon  release  for  each  replication 

Index  of  submarine  attacks  on  main  body  ships 


3. 14159 


Index  used  in  subroutine  FIND  to  isolate  sections  of 


the  subroutine  after  the  first  call  on  the  subroutine 
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Name 


Definition 


SONAR(I,  J) 


SONAV(I) 


SSV(I,  J) 
SSV(I,  1) 
SSV(I,  2) 
SSV(I,3) 
SSV(I,  4) 
SSV(I,  5) 
SSV(I,  6) 
SSV(I,  7) 

SUBKIL 

SV(I,  J) 


SV(I,  1) 

SV(I,  2) 
SV(I,  3) 
SV(I,  4) 
SV(I,5) 
SV(I,6) 
SV(I,7) 
SV(I,  8) 


The  generated  detection  range  ofrthe  J1     screen  ship 
for  the  I"1  replication  in  a.  run 


The  average  of  generated  detection  ranges  of  the  I 


th 


screen  ship  in  a  run 


rth 


Information  vector  of  the  I       submarine  in  a  replication 
Number  of  the  I      submarine 


Submarine  course  (in  radians) 


Submarine  speed 


Submarine  X  location 


Submarine  Y  location 


Time  X_Y  location  fixed 


Number  of  times  detected  by  screen  ships 

Number  of  times  the  penetrating  submarines  detected 

Information  vector  of  the  I       screen  ship  in  a  repli- 


cation 


th 


Number  of  the  I       screen  ship 

Screen  ship  course  (in  radians) 

Screen  ship  speed 

Screen  ship  X  Location 

Screen  ship  Y  location 

Screen  ship  sonar  type 

Screen  ship  generated  detection  range 


Patrol  area  radius 
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Name 

SV(I,  9) 
SV(I,  10) 
SV(I,  11) 
SV(I,  12) 
SV(I,  13) 

TAR(I,  J) 
TAR(I,  1) 
TAR(I,  2) 
TAR(I,  3) 
TAR(I,  4) 
TAR(I,  5) 

TEMP(I) 


Definition 


X  location  of  center  of  assigned  patrol  area 
Y  location  of  center  of  assigned  patrol  area. 
Time  of  next  course  change 
Number  of  submarines  detected  in  replication 


Time  X-Y  location  fixed 


rth 


Information  vector  for  the  I      main  body  ship 

Number  of  the  I**1  main  body  ship 

Main  body  ship  course  (in  radians) 

Main  body  ship  speed 

Main  body  ship  X  location 

Main  body  ship  Y  location 

Used  to  aid  in  ordering  event  calendar  as  a  temporary 

storage  location  for   an  event 
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by  some  screening  ships,    and  some  penetrating  submarines.      The  model 
is  designed  for  use  as  an  aid  in  improving  the  ability  of  a  proposed  screening 
tactic  in  the  detection  of  a  penetrating  submarine.     A  systematic  procedure 
to  improve  a  screen's   effectiveness  against  a  known  submarine  threat  is 
demonstrated,    and  an  example  problem  is  worked  using  this  procedure. 
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